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Abstract. In minimal supergravity (mSugra), the neutrino sector is related to the slepton 
sector by means of the renormalization group equations. This opens a door to indirectly test 
the neutrino sector via measurements at the LHC. Concretely, for the simplest seesaw type-I, 
we present the correlations between seesaw parameters and ratio of stau lepton flavour violating 
(LFV) branching ratios. We find some simple, extreme scenarios for the unknown right-handed 
parameters, where ratios of LFV rates correlate with neutrino oscillation parameters. On the 
other hand, we scan the mSugra parameter space, for both seesaw type-I and II, to find regions 
where LFV stau decays can be maximized, while respecting low-energy experimental bounds. 
We estimate the expected number of events at the LHC for a sample luminosity ol C = lOOfb"'^. 



1. Introduction 

Many theoretical models aim at explaining what neutrino experiments [1-5] have demonstrated: 
that neutrinos have small mass. Among the proposals to explain current neutrino data [6,7], 
the most popular ones are based in the seesaw mechanism [8-13] . Although this kind of models 
are not directly testable, they can be indirectly tested if we assume a supersymmetry (SUSY) 
framework with universal boundary conditions (like mSugra), as the Renormalization Group 
Equations (RGE's) relate the neutrino Yukawa couplings to the soft SUSY breaking slepton 
mass parameters. Section 2 shows the correlations between stau LFV branching ratios (BR's) 
and neutrino parameters. In the case of type-I seesaw, we only present the correlations obtained 
assuming degenerate heavy neutrinos. For a more complete analysis, see reference [14]. Similar 
correlations as the ones obtained in the simplest type-I SUSY seesaw can be obtained in a 
SU{5) inspired type-II SUSY seesaw model [15]. For more details about the correlations, see 
reference [16]. Section 3 presents an analysis over the mSugra parameter space of the stau 
LFV decays and the estimated maximum number of events for the opposite-sign dilepton signal 
X2 ~^ Xi f^'^j for both type-I and II seesaw mechanisms. 



2. Correlations for type-I seesaw 

In section 2.1 we will obtain analytical expressions for the correlations between neutrino 
parameters and stau LFV BR's, while in section 2.2 we will check the validity of the analytical 
estimates via a full numerical calculation. 



2.1. Analytical correlations 

We briefly introduce the simplest type-I SUSY seesaw mechanism. Its particle content is the 
same as in the MSSM, but enlarged by three right-handed neutrino superfields Nf. The leptonic 
part of the superpotential is 

W = Y^%H^'^ + Yl%HuN^ + MiNfNf , (1) 

where and Y^, denote the charged lepton and neutrino Yukawa couplings, while Mj are the 
Majorana mass terms (of unspecified origin) of the "right-handed" neutrino superfields. From 
now on, we will work in the basis where both the Majorana mass matrix of the right-handed 
neutrinos and the charged lepton Yukawa coupling matrix are diagonal: M/j and 1^, respectively. 

The effective mass matrix of the mostly left-handed neutrinos is given in the usual seesaw 
approximation as 

m. = -^yJ-M^i-y.. (2) 
This complex symmetrix matrix is diagonalized by the leptonic mixing matrix U [12] 

rhi, = ■ nil, ■ U , (3) 
where rhi, = diag(mi). We can parametrize the neutrino Yukawa matrix as [17] 

Y^ = V2—\/^R^/ih^U\ 
Vu 

where rhi, and Mr are diagonal matrices with the light neutrino mass eigenvalues rui and the 
heavy neutrino mass eigenvalues Mj, respectively; U is the leptonic mixing matrix and R is a 
complex orthogonal matrix. 

The LFV branching ratios of the charged sleptons are related to the off-diagonal elements of 
their mass matrix. Neglecting left-right mixing in the charged slepton sector, these correspond 
to the off-diagonal elements of the slepton soft SUSY breaking mass parameter. In the general 
MSSM, the soft SUSY breaking parameters are free. However, if we assume an mSugra-like 
framework as the origin of SUSY breaking, soft SUSY breaking masses are diagonal at the GUT 
scale and they only get non-zero off-diagonal elements at the electroweak (EW) scale by means 
of the RGB's. If the mixing between different flavour eigenstates of the charged sleptons is 
small (what we call the small-angle approximation), the LFV BR's for the decays of the charged 
sleptons can be approximated in terms of the off-diagonal soft breaking mass parameters 
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(5) 

In the leading-log approximated solutions to the RGE's, the left-slepton LFV decays are 
proportional to 

(6) 

where L is deflned as 

and Mx is the GUT scale. Note that, in the leading-log approximation, the off-diagonal elements 
of the right-slepton soft breaking mass are zero, so that we will restrict our analysis to "left- 
handed" sleptons. Using the parametrization defined in equation (4), the left-slepton LFV 
decays can be correlated to neutrino parameters 
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Taking ratios of LFV bracking ratios, we can eliminate the dependence on the SUSY parameters, 
so that we can establish clean correlations between "left-handed" slepton LFV BR's and neutrino 
parameters only. This way, for example, we can estimate that the ratio of the heaviest stau LFV 
decays is related to neutrino parameters as 
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We can use equation (9) to obtain precise correlations between observables for different neutrino 
scenarios. If we assume that heavy neutrinos are degenerate, R is real and that U fulfills the 
tribimaximal (TBM) ansatz [18], then 



^13 
'23 



2()-;?2 - nil) 
|3m3 — 2m2 — mi\ 



(10) 



Table 1 shows the form of equation (10) and its numerical values, for different neutrino scenarios. 
For a general analysis of different neutrino hierarchies, dependence on the overall light neutrino 



Table 1. Parameter for the case of degenerate heavy neutrinos, assuming R being real and 
TBM mixing for light neutrinos. Each column label corresponds to a different light neutrino 
scenario: strict normal hierarchical (SNH), strict inverse hierarchical (SIH), quasidegenerete 
normal hierarchical (QDNH) and quasidegenerete inverse hierarchical (QDIH) light neutrinos. 
Row label ANALYTICAL shows the analytical form of the parameter r2^; row label BFP shows 
its value when fixing neutrino mass splittings to their best fit point value [6, 7]; and row label 3a 
shows its value when considering the 3a allowed range for the neutrino mass-squared differences. 
Note that a = Amg/\Am1\ is the ratio of the solar over the atmospheric mass splitting and 
(Ta = Am2/|Am2| is the sign of the atmospheric mass splitting. 



SNH 



SIH 



ANALYTICAL 
BFP 

3(7 




1.7 X 10-2 



^13 _ 
'23 ~ 

(41? 



1.1 X 10" 



(r^f)2 G [0.92, 3.7] x lO'^ {rgf G [0.48, 3.3] x 10" 



QDNH 



QDIH 



ANALYTICAL 
BFP 

3(7 



2a 



j,13 ^ 

23 3(7\-\-a 

(rM)2 G [1.8, 12] X 



4.4 X 10" 



10" 



„13 
'23 

(ri|)2 = 4.6 X 10-"^ 
(ri)2 G [1.9, 13] X 10" 



2a 



3(Ta 



mass scale, departure from tribimaximality, both for degenerate and strongly hierarchical heavy 
neutrinos, see reference [14]. 

2.2. Numerical correlations 

In order to check the analytical correlations obtained in section 2.1, we have performed a 
numerical calculation using the program package SPheno [19]. Although we have checked 
that the validity of the analytical correlations is independent of the mSugra parameters, we 
only present the results for some specific benchmark points. A more datailed study over the 
mSugra parameter space is presented in section 3. The numerical procedure followed to fit 



neutrino data consists in first fixing the neutrino Yukawa couplings at the GUT scale by means 
of equation (4) and assuming certain neutrino scenarios as input parameters. The RGE's are 
numerically run with SPheno, so that the light neutrino masses are calculated at the EW scale 
as output parameters. An iterative process varies the input parameters at the GUT scale in 
order to obtain the desired output parameters at the EW scale. 

Figure 1 shows the dependence of LEV processes (low energy LEV processes in the left panel 
and stau LEV decays in the right panel) on the common mass Mr of the degenerate heavy 
neutrinos, for the mSugra benchmark point SPS3 [20]. We have checked that our analytical 
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Figure 1. Branching ratios for li Ij j and Zj — > 31 j (left), and f2 ex? and f2 /ix? 
(right) for the standard point SPS3 versus Mji, assuming degenerate right-handed neutrinos. 
Neutrino oscillation parameters have been fixed to the best fit values for the mass splittings [6, 7], 
with exact TBM neutrino angles. The absolute neutrino mass scale has been fixed to mi = 0. 
The coloured region in the right-side plot is excluded from the current experimental limit on 
BR(Ai ^ 67). 

estimate of ratios of stau LEV branching ratios is in very good agreement with the numerically 
calculated ratio, as long as the magnitude of the stau LEV BR's are not larger than 10% 
(where the small-angle approximation is no longer valid). It is important to note that current 
experimental upper bounds on low energy LEV processes put a severe constraint in the maximum 
possible value of stau LEV decays. Given the current experimental limit on BR(// — > 67) [21] 

BR{p 67) < 1.2 • 10"^^ , (11) 

the maximum possible value for BR(f2 MX?) is a few percent for SPSS. 

On the contrary, for benchmark point SPSla' [22], the allowed maximum value for BR(f2 — > 
HXi) is only 10^^. The reason for this difference between points SPS3 and SPSla' is that the 
SUSY mass spectrum is nearly two times heavier in the former. As BR(/x — > 67) depends 
inversely on the SUSY spectrum [23], 

BR(^^e7) oc (12) 
"^SUSY 

its value is several hundred times smaller for SPS3, thus providing a weaker constraint on f2 
LEV BR's. 

The dependence of BR(;U — > 67) on the SUSY spectrum is shown in figure 2. Here we plot 
BR(// 67), BR(-f2 — > fJ-Xi) and BR(f2 6X?) versus the mass of f2. The range of nir^ has 
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Figure 2. LFV BR's obtained as function of the heavier scalar tau mass. The fuh hne (red) 
is 10^^ X BR(/x ^67), the dashed hne (blue) BR(f2 fJ-Xi) and the dot-dashed line (green) 
is BR(f2 — > exi)- Data obtained for SPSla with parameters varied along the corresponding 
"slope". Right-handed neutrino mass is fixed to Mr = 3 x 10^^ GeV. Neutrino oscillation 
parameters have been fixed to the best fit values for the mass splittings [6, 7], with exact TBM 
neutrino angles, and the absolute neutrino mass scale has been fixed to nii = 0. The black line 
is the current upper limit on BR(^ — > 67). 



been obtained by varying the SPSla parameters along the "slope" [20]. Note that BR(/x 67) 
decreases as the SUSY spectrum gets heavier, dropping below its current experimental bound for 
TUf^ larger than about 250 GeV. In contrast, f2 LFV BR's increase for increasing rrif^, as charged 
left-sleptons become increasingly degenerate and thus their mixing increases. Nevertheless, note 
also that the ratio of f2 LVF BR's remains constant, in agreement with the analytical estimate, 
as long as the magnitude of each individual BR is not too large (and hence, the small-angle 
approximation is valid). 

3. Scan 

In section 2 it has been shown that it is possible to indirectly test the seesaw mechanism at 
the LHC by means of correlations between the ratio of f2 LFV BR's and neutrino parameters. 
Although the magnitudes of f2 LFV BR's have been studied for two specific SUSY benchmark 
points (SPSla' and SPSS), a more general study over the mSugra parameter space is necessary. 
Thus, we can identify regions of the mSugra parameter space that maximize the magnitude of 
f2 LFV BR's. Besides this, we have estimated the maximum number of events of the opposite- 
sign dilepton signal X2 ~^ XifJ-Ti which can be searched for at the LHC. Note that a complete 
Monte Carlo analysis would be needed, but this is out of the scope of this work. The numerical 
procedure consists in scanning over the mQ-mi/2 plane, for fixed values of other mSUGRA 
parameters, using the program package SPheno [19]. For each point in this plane, we perform 
a maximization of BR(^ — > 67), so that it is as close as possible to its current experimental 
upper bound, given in equation (11). To fit neutrino data we have followed the same iterative 
procedure described in section 2.2. In our analysis we have always fitted light neutrino mass 
splittings to their best fit point values [6,7], under a strictly normal hierarchy (m-i = 0) and 
their mixing to be TBM. To calculate the number of events of the opposite-sign dilepton signal 
X2 ~^ Xi A*''") we have estimated the total production cross section of X2 at leading order with 
the package Prospino [24-28]. 



3.1. Type-I seesaw 

In order to simplify our numerical analysis in type-I seesaw, we have only considered a specific 
neutrino scenario, assuming that right-handed neutrinos are degenerate and the matrix R is 
real. Figure 3 shows the magnitude of BR(f2 A*Xi) (l^ft column) and BR(f2 ex?) (right 
column) in the mo — plane, for three combinations of mSugra parameters: the first row 
corresponds to what we call our standard mSugra point, defined by 

tan/? =10, ^0 = 0, /x>0; (13) 

the second row corresponds to the same standard mSugra point, but with tan/3 = 30; the third 
row corresponds to the same standard mSugra point, but with = —300 GeV. One important 
feature is tliat there is a region in the mo — mi/2 plane (generally, small mo and large mi/2 
values) where f2 LFV BR's can be as large as 10%. This is the parameter space region where 
the constraint imposed by BR(/i ^ 67) is supressed due to a strong cancellation between its 
neutralino-charged slepton and the chargino-sneutrino amplitude contributions. Note that, as 
previously shown in section 2.1, the BR(f2 ex?) is normally smaller than the BR(f2 /^Xi) 
by a factor 1.7 x 10~^. While f2 LFV BR's are decreased by larger values of tan/3, they are 
practically not affected by ^o- We have checked that the f2 LFV BR's are not sensitive to the 
sign of neither Aq nor /x. 

In figure 4, we have plotted the production cross section cr(x2) at leading order times the 
BR of X2 go™g to the opposite-sign dilepton signal Xiywr as a function of mi/2) for different 
values of mo. We have fixed the rest of the mSugra parameters to our standard point defined 
in equation (13). For mo ^ 100 GcV and mi/2 ~ [450, 600] GeV and assuming a luminosity 
C = WOfb^^, the number of events of the opposite-sign dilepton signal X2 Xi A*''" can be of 
the order of 10^. 

3.2. Type-II seesaw 

In order to simplify our numerical analysis in type-II seesaw, we have considered only the case 
in which Ai = 2 x 10~^ and A2 = 0.5. We have explicitly checked that Ai plays no role in 
LFV. On the other hand, larger values of A2 generally increase the LFV signal, as BR(// — e 7) 
decreases for larger A2. Figure 5 shows the same as figure 3, but for type-II seesaw. Note that 
the region where f2 LFV BR's is smaller and a deformed with respect the one for type-I seesaw. 
The reason for this is that the addition of non-gauge-singlet states in type-II seesaw increases 
the dependence on the renormalization scale of the neutrino Yukawa coupling and also affects 
the SUSY spectrum (and thus the region where BR(;U ^67) is strongly suppressed). 

In figure 6 it is plotted the same as in figure 4, but for type-II seesaw. Assuming a luminosity 
oi £. = 100/6~^, there can be a maximum number of events of the order of 10^ for mo ~ 100 
GeV and mi/2 ~ [600, 800] GeV. 

4. Conclusions 

We have shown that neutrino parameters can be indirectly tested at the LHC by measuring the 
ratio of f2 LFV BR's in type-I and II SUSY seesaw models if mSugra is assumed. We have 
performed a numerical analysis of the absolute values of f2 LFV BR's and we have estimated 
the maximum number of events that can occur at the LHC. We have shown that there exist 
regions of the mSugra parameter space where the number of events of the opposite-sign dilepton 
signal X2 Xi A*''" can be as much as of the order of 10^. 
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Figure 3. BR(f2 — > ^x?) (l^ft column) and BR(f2 — > ex?) (right column) in the mo — m.x/2 
plane, for the type-I seesaw. We impose BR(^ ^ 6 7) < 1.2 • 10-i\ and we consider different 
choice of mSugra parameters: /i > 0, tan/3 = 10 and = GeV (first row); > 0, tan/3 = 30 
and = GeV (second row); fJ- > 0, tan/3 = 10 and Aq = —300 GeV (third row). 
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Figure 4. Production cross section (at leading order) of X2 times BR(x2 — ^ Xi f^-'^) versus mi/2 
for mo = 100 GeV (red), 200 GeV (green), 300 GeV (blue) and 500 GeV (yellow). We take 
a standard choice of parameters: ^ > 0, tan /? = 10 and ^Iq = GeV, for type-I seesaw, for 
BR(/i ei) < 1.2 • 10"^^ 
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Figure 5. BR(f2 fJ-Xi) (l^ft column) and BR(f2 — > ex?) (right column) in the mo — m.x/2 
plane, for the type-II seesa'w. We impose BR(^ ^ 67) < 1.2 • 10-i\ and "we consider different 
choice of mSugra parameters: /i > 0, tan/3 = 10 and = GeV (first row); > 0, tan/3 = 30 
and Ao = GeV (second ro"w); > 0, tan/3 = 10 and Aq = —300 GeV. (third ro"w) 
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Figure 6. Production cross section (at leading order) of X2 times BR(x2 — ^ Xi f^-'^) versus mi/2 
for mo = 100 GeV (red), 200 GeV (green), 300 GeV (blue) and 500 GeV (yellow). We take 
a standard choice of parameters: /u > 0, tan/3 = 10 and Aq = GeV, for type-II seesaw, for 
BR(/i ^ e + 7) < 1.2 • 10"^^ 



SUSY 02, DESY Hamburg, June 17-23, 2002 vol 1 ed Nath P, Zerwas P M and Grosche C (Hamburg: 
DESY) pp 217-26 [Preprint hep-ph/0211145) 
[28] Plehn T 2005 Czech. J. Phys. 55 B213-B220 {Preprint hep-ph/0410063) 



